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ABSTRACT 



We report the discovery of the low-density, transiting giant planet WASP-31b. The planet is 0.48 Jupiter masses and 1.55 Jupiter radii. 
It is in a 3.4-day orbit around a metal-poor, late-F-type, V = 1 1.7 dwarf star, which is a member of a common proper motion pair. In 
terms of its low density, WASP-31b is second only to WASP- 17b, which is a more highly irradiated planet of similar mass. 

Key words, binaries: eclipsing - planetary systems - stars: individual: WASP-31 



1. Introduction 

To date, 107 transiting extrasolar planets have been disco verecQ, 
the majority of which are gas giants in short orbits. The radii of 
a subset of these exoplanets are larger tha n predicted by stan- 
dard models of irra diated gas giants (e.g., Burrows et al. 2007; 
i Fortnev et alJ2007l) . including TrES-4b dMandushev et aLl2 007; 
Soz zetti et aU2009h WASP- 12b dHebb et alj|2009h . and WASP- 
17b (lAnderson et alJ[2Q10L 1201 lbT A number of mechanisms 
hav e been proposed as po tential solutions to the radius anomaly 
(see IFortnev et alJ d2010l) for a review), each of which involves 
either injecting heat into the planet from an external source or 
slowing heat loss from the planet. 

One such mechanism is the dissipation of energy within a 
planet as heat during the tidal circularisation of an eccentric orbit 
dBodenheimer et al.ll200il: iGu et alJl200l iJackson et al]l2008t 
Ibgu i & Burrowsl l2009). Such studies suggest that tidal heating 
may be sufficient to explain the large radii of even the most 
bloated exoplanets, though we would have to be observing some 
sy stems at very special times. A high heating rate, as suggested 
bv iLeconte et alJ d2010h . would mean most tidal energy is radi- 
ated away by the age typical of the very most bloated planets 
(a few Gyr) and so could not have played a significant role in 



* Based in part on observations made with the HARPS spectro- 
graph on the 3.6-m ESO telescope (proposal 085.C-0393) and with the 
CORALIE spectrograph and the Euler camera on the 1.2-m Euler Swiss 
telescope, both at the ESO La Silla Observatory, Chile. 
** The photometric time-series and radial-velocity data used in 
this work are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsarc.u-strasbg.fr/viz-bin/qcat7J/A-l-A/531/A60 

1 2010 Nov 25, http://exoplanet.eu 



their observed bloating. However, the current uncertai nty in tidal 
theory allows for a wide range of heating rates (e.g. Ibgu i et all 
201 1). Though most studies have consi dered a transient p hase of 
tidal heating, ongoing tidal heating (e.g. Ibg ui et al.l2010l) would 
occur if an additional companion conti nues to excite the orbital 
eccentricity of the bloated planet (e.g.[Mardling 2 0071) . 

iBurrows et al.l d2007l) proposed that enhanced opacities 
would retard the loss of internal heat and thus slow contraction 
of bloated planets. They suggested that enhanced opacities may 
arise due to the strong optical and UV irradiation of short-orbit, 
gas giants that could alter their atmospheres, producing thick 
hazes, absorbing clouds and non-equilibrium chemical species 
(e.g. tholins or polyacetylenes). 

The bloated planets are all very strongly irradiated by their 
host stars, and a small fraction of stellar insolation energy would 
be sufficient to account fo r the observed degrees of bloating. 
iGuillot & Showman] d2002l) suggested that the kinetic energy 
of strong winds, induced in the atmosphere by the large day- 
night temperature contrasts that result from tidal locking, may 
be transported downward and deposited as thermal energy in the 
deep interior. However, a mechanism to convert the kinetic en- 
ergy i nto ther mal energy would s t ill be r equired. iLi & Good man 
(l2010h and lYoudin & Mitchelll d2010h found that turbulence 
is efficient at dissipati ng kinetic energy . Magnetic drag on 
weakly ionized winds dPerna et al.l 12010b and Ohmic heating 
dBatvgin & Stevenson! 1201 Oh are alternative mechanisms. The 
non-bloated planets are also highly irradiated. Hence, such a 
mechanism would either have to act more efficiently on the 
bloated planets, or some other property must counteract its ef- 
fect. O ne such possibility is th e presence of a mas sive core. 
Indeed, IGuillot et all d2006l) and IBurrows et ail d2007l) found a 
correlation between the core masses required to reproduce the 
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Fig. 1. WASP-South discovery light curve. Upper panel: 
Photometry folded on the orbital period of P — 3.4 d. Points 
with error above three times the median error (0.012 mag) were 
clipped for display purposes. Lower panel: Photometry folded 
on the orbital period and binned in phase (A0 = 0.025), with the 
transit model generated from the parameters of Table [4] superim- 
posed. 

Table 1. Proper motions of WASP-31 and its visual companion 



Star 


Catalogue 


yU RA (mas) 


//Dec (mas) 


WASP-31 


UCAC3 


-28.2 ± 1.3 


-0.4 ± 1.6 


WASP-31 


PPXML 


-25.0 ± 2.3 


-0.1 ± 2.4 


Companion 


UCAC3 


-33.1 ±3.1 


+ 1.0 ± 3.8 


Companion 


PPXML 


-28.5 ± 4.2 


+ 1.6 ±4.2 



observed radii of known exoplanets and the metallicities of their 
host stars. 

In this paper, we present the discovery of the bloated, tran- 
siting, giant planet WASP-31b. Compared to the ensemble of 
known short-period planets, WASP-31b is moderately irradiated 
by its low-metallicity host star. 

2. Observations 

WASP-31 is a V = 11.7, F7-8 dwarf star located in the con- 
stellation Crater. WA SP-31 has been observed by WASP-South 
(Pollacco et al. 2006) during the first five months of each year 
since t he start of full-scale operati ons (2006 May 4). A transit 
search (ICollier Cameron et al.ll2006h of the resulting 24614 us- 
able photometric measurements (Figure [TJ found a strong, 3.4-d 
periodicity. 

WASP-31 is a visual double with a V ~ 15.8 star (2MASS 
11174477-1903521) approximately 35" away. The 2MASS 
colours of the companion suggest that it is a mid-to-late K- 
type star. The proper moti ons for the t wo stars listed in the 
PPMXL dRoeseretalj2010l) and UCAC3 dZacharias et alj2 010) 
catalogues suggest that this is a common proper motion pair 
(Table [TJ. The companion is blended with WASP-31 in the 
WASP images, so we corrected the WASP photometry for this 
contamination to prevent dilution of the transit. 

Using the CORALIE s pectrograph mounte d on the 1.2- 
m Eu ler-Swiss telescope (Bara nne et al.l Il996t lOueloz et alj 
l2000bl) . we obtained 34 spectra of WASP-31 during 2009 and a 



further 13 spectra during 2010. As pressure variations can cause 
CORALIE to drift on short timescales, we calibrated the stel- 
lar spectra by obtaining simultaneous spectra of a thorium-argon 
lamp. In April 2010, we obtained an additional 10 spectra with 
the HARPS spectrograph mounted on the 3.6-m ESO telescope. 
As HARPS is stable at the 1ms" 1 night -1 level, we obtained a 
calibration at the start of each night of observations. This avoids 
contamination of the stellar spectra by the thorium-argon lamp. 
The CORALIE measurement taken at BJD = 2454971.548671 
and the HARPS measurement taken at BJD = 2455 299.716991 
were both affected by cloud cover. 

The typical signal-to-noise ratio (S/N) per pixel at 550 nm 
is 18 for the CORALIE spectra, with exposure times of 30 min, 
and 26 for the HARPS spectra, with exposure times of 15 min. 
Radial-velocity ( RV) measurements were computed by w eighted 
cross-correlation dBaranne et alJ U996: Pep e et aUl2005h with a 
numerical G2-spectral template. RV variations were detected 
with the same period found from the WASP photometry and with 
a semi-amplitude of 58 m s , consistent with a planetary-mass 
companion. The RV measurements are plotted in the upper panel 
of Figure [2] 

To test the hypothesis that the RV variations are due to spec- 
tral line distortions caused by the presence of cool s tellar spots or 
a blen ded eclipsing binary, a line-bisector analysis dOueloz et al.l 
I200TJ: ISantos et alj|2002l) of the CORALIE and HARPS cross- 
correlation functions was performed. The lack of correlation be- 
tween bisector span and RV (Figure [2 lower panel), especially 
for the high-precision HARPS measurements, supports the iden- 
tification of the transiting body as a planet. As an additional 
test for spots or a blend, we also computed CORALIE RVs by 
weighted cross-correleation with a numerical K5-spectral tem- 
plate. The amplitude and phase of the RV variations are the 
same within errors, irrespective of the choice of cross-correlation 
mask, as would be expec ted if the variations w ere caused by the 
presence of a planet (e.g. lHuelamo et ai1l2008l) . 

To refine the system paramters, we obtained high-S/N tran- 
sit photometry. Photometric follow-up observations of WASP-3 1 
were obtained with the LCOGlQ 2.0-m Faulkes Telescope North 
(FTN) on Mt. Haleakala, Maui on the night of 2010 Feb 26. The 
fs03 Spectral Instruments camera was used with a 2 x 2 bin- 
ning mode, giving a field of view of 10' x 10' and a pixel scale 
of 0.303" pixel -1 . The data were taken through a Pan-STARRS 
z filter, with the telescope defocussed to minimise flat-fielding 
errors and to allow 60-s exposure times to be used without satu- 
rating. 

The data were p re-processed using the WASP Pipeline 
dPollacco et al.l [20060 to perform masterbias and flat construc- 
tion, debiassing and flatfielding. Due to the very low dark cur- 
rent of the fs03 Fairchild CCD (< 0.0001 e" pix" 1 s" 1 ), dark sub- 
traction was not performed. Aperture photometry was performed 
using DAOPHOT within the IRAF environment using an aper- 
ture with a radius of 1 1 pixels. Differential photometry was then 
performed relative to 20 comparison stars that were within the 
FTN field of view (Figure [3]). The RMS of the residuals about 
the best-fitting model (Section|4| was 1.4 mmag. 

On 2010 April 15 we obtained 4.1 hours of photometry in 
the Gunn r filter with the CCD camera on the Euler-Swiss tele- 
scope, covering from 40 min before the start of transit until 55 
min after it ended. The conditions were variable, with seeing of 
0.6-1.7" and an airmass range of 1.15-1.34. Euler now employs 
absolute tracking to keep the stars on the same pixels during a 
whole transit. By identifying point sources in each image and 



2 http://lcogt.net 
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Fig. 2. Upper panel: Spectroscopic orbit of WASP-31, as il- 
lustrated by radial velocities from CORALIE (blue circles) and 
HARPS (red squares). The best-fitting Keplerian model, gener- 
ated from the parameters of Table |4] is overplotted as a solid 
line. An RV taken at BJD = 2 455 168.8468, depicted in the plot 
by an open circle, fell during tra nsit. As we did no t treat the 
Rossiter-McLaughlin effect (e.g.. lOueloz et al]|2000al) . we ex- 
cluded this measurement from our combined analysis. Lower 
panel: A lack of correlation between bisector spans and radial 
velocities rules out a blended eclipsing binary or starspots as 
the cause of the photometric and spectroscopic variaitions. We 
adopted uncertainties on the bisector spans twice the size of 
those on the radial velocities. For both plots, the centre-of-mass 
velocity, y = -124.92 ms _I , was subtracted from the radial ve- 
locities and the Keplerian model. 



matching them with a catalogue, the image centre is calculated. 
Drifts from the nominal position are then corrected by adjusting 
the telescope pointing between exposures. 

After bias-subtracting and fiat-fielding the images, we per- 
formed aperture photometry. The flux was extracted for all stars 
in the field and the final light curve (Figure 0) was obtained by 
differential photometry of the target and a reference source ob- 
tained by combining the 4 brightest reference stars. The RMS 
of the residuals about the best-fitting model (Section|Ui was 2.6 
mmag, which was limited by the number of available reference 
stars. 
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Fig. 3. High-S/N transit light curves. The upper observations 
(blue circles) were obtained by FTN, using a Pan-STARRS z fil- 
ter, on 2010 Feb 26. The lower observations (red diamonds), off- 
set in relative flux by 0.026 for display, were obtained by Euler, 
using a Gunn r filter, on 2010 Apr 15. The best-fitting transit 
models generated from the parameters of Table [4] are overplot- 
ted. 



3. Stellar parameters 

The individual HARPS spectra of WASP-31 were co-added 
to produce a single spectrum with an average S/N of around 
100:1. The analysi s was performed using the methods given in 
iGillon et alj d2009t) . The H a line was used to determine the ef- 
fective temperature (T e ff), while the Na i D and Mg i b lines 
were used as surface gravity (logg*) diagnostics. The parame- 
ters obtained from the analysis are listed in Table[2] The elemen- 
tal abundances were determined from equivalent width measure- 
ments of several clean and unblended lines. A value for micro- 
turbulence (&) was determined from Fe i using the method of 
Magain ( 1984). The quoted error estimates include that given by 
the uncertainties in T e g, logg, and tf t , as well as the scatter due 
to measurement and atomic data uncertainties. 

The sky-projected stellar rotation velocity (v sin 7) was de- 
termined by fitting the profiles of several unblended Fe i lines. 
We assumed a value for macrotur bulence (v ma c) of 5.2 ± 0.3 km 
s _1 , based on the tabulation bv iGravl d2008h . and we used an 
instrumental FWHM of 0.06 ± 0.01 A, determined from the tel- 
luric lines around 6300A. A best-fitting value of v sin / - 7.6 + 
0.4 km s -1 was obtained. However, recent work bv feruntt et all 
(2010) suggests a lower value for macroturbulence of v mac = 4.2 
+ 0.3 km s~ 1 which yields a slightly higher v sin / = 8 . 1 +0.4 km 
s . We therefore adopt the average of these two determinations, 
v sin / = 7.9 + 0.6 km s _1 , with the uncertainty being the quadra- 
ture addition of the individual uncertainties. If v mac = km s _1 , 
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Table 2. Stellar parameters from the spectroscopic analysis 



Parameter 


Value 




6300 ± 100 K 


logg. 


4.4 + 0.1 (cgs) 




1.4 ±0.1 km s _I 


v sin/ 


7.9 ±0.6 km s _1 


[Fe/H] 


-0.20 ± 0.09 


[Na/H] 


-0.24 ± 0.04 


[Mg/H] 


-0.11 ±0.06 


[Si/H] 


-0.13 ±0.07 


[Ca/H] 


-0.03 ± 0.08 


[Sc/H] 


-0.07 ± 0.06 


[Ti/H] 


-0.10 ±0.09 


[V/H] 


-0.16 ±0.09 


[Cr/H] 


-0.20 ± 0.08 


[Mn/H] 


-0.45 ±0.11 


[Ni/H] 


-0.25 ± 0.08 


log A(Li)[LTE] 


2.82 ± 0.08 


log A(Li)[NLTE] 


2.75 ± 0.08 


M, 


1.15 ± 0.08 M Q 


R, 


1.12 ± 0.15 R Q 


R.A. (J2000) 
Dec. (J2000) 

Knag 
Jmag 
^mag 
^rriag 

USNO-Bl.O 
2MASS 


ll h 17 m 45.35 s 
-19°03'17.3" 

11.7 ±0.2 
10.91 ± 0.02 
10.71 ± 0.02 
10.65 ± 0.03 
0709-0239208 
11174536-1903171 



Note: NLTE Lithium value 
Mass and Radius estimate 



using correction of Carlsson et al.l d 1 994T) . 
using the lTorres et al] f2010h calibration. 



Table 3. Limb-darkening coefficients 



Light curve (band) 


a\ 


a 2 


«3 


a 4 


WASP/Euler (RC) 


0.430 


0.488 


-0.254 


0.020 


FTN (Pan-STARRS z) 


0.520 


0.060 


0.100 


-0.107 



then a value of v sin / = 8.7 + 0.4 km s 1 is found, which is the 
upper-limit of the sky-projected rotation velocity. 



4. Combined analysis 

The WASP, FTN and Euler photometry were combined 
with the CORALIE and HARPS radial velocities in a 
simul taneous Markov-chain Mo n te Carlo (MCMC) analy- 
sis (ICollier Cameron et al.l 120071: iPollacco et all 12008b . The 
transit light curve w as modeled using the formulation of 
iMandel & Agoll (|2002|) with the assumption that R pi <s:R t . Limb- 
darkening was accounted for using a four-coefficient nonlinear 
limb-darkening model, using fixed coefficients (Table [3]) appro- 
priate to the passbands and interpolated in effective temper ature, 
surfac e gravity and metallicity from the tabulations of lOaretl 
(I2000h . 

The transit light curve is parameterized by the epoch of 
midtransit 7*o, the orbital period P, the planet-to-star area ratio 
(R p \/R*) 2 , the approximate duration of the transit from initial to 
final contact Tu, and the impact parameter b — a cos i/R* (the 
distance, in fractional stellar radii, of the transit chord from the 
star's centre). The radial-velocity orbit is parameterized by the 
stellar reflex velocity semi-am plitude K*, the systemi c velocity 
y, and y[e cos u> and ye sin w (lAnderson et al.11201 lab , where e 
is orbital eccentricity and co is the argument of periastron. 



The linear scale of the system depends on the orbital sepa- 
ration a which, through Kepler's third law, depends on the stel- 
lar mass M,. At each step in the Markov chain, the latest val- 
ues of stellar density p*, effective temperature r e (j and metal- 
licity fFe/HI are in put in to the empirical mass calibration of 
lEnoch et alT( 20101) to obtain M*. The sh apes of the transit light 
curves dSeager & Mallen-Ornelasll2003l) and the radial-velocity 
curve constrain p„, which combines with M, to give R„. T e s and 
[Fe/H] are proposal parameters constrained by Gaussian priors 
with mean values and variances derived directly from the stellar 
spectra (Table |2]l. 

As the planet-star area ratio is constrained by the measured 
transit depth, R p \ follows from R». The planet mass M v \ is cal- 
culated from the measured value of K\ and M*; the planetary 
density p p \ and surface gravity log g p \ then follow. We also calcu- 
late the blackbody equilibrium temperature r eq i (assuming zero 
albedo and efficient redistribution of heat from the planet's pre- 
sumed permanent day side to its night side), the transit ingress 
and egress durations, T\2 and T34, and the orbital semi-major 
axis a. 

At each step in the MCMC procedure, model transit light 
curves and radial-velocity curves are computed from the pro- 
posal parameter values, which are perturbed from the previous 
values by a small, random amount. The x 2 statistic is used to 
judge the goodness of fit of these models to the data and a step is 
accepted if x 2 is lower than for the previous step. A step with 
higher x 2 is accepted with a probability exp(-A^- 2 /2). In this 
way, the parameter space around the optimum solution is thor- 
oughly explored. To give proper weighting to each transit and 
RV data set, the uncertainties are scaled at the start of the MCMC 
so as to obtain a reduced x 2 of unity. We allow for a system- 
atic instrumental offset, AyeARPS, between the CORALIE and 
HARPS spectrographs. 

From an initial MCMC fit for an eccentric orbit, we found 
e = 0.02 7 + ggl with a 3-cr upper limit of 0.13. The F-test ap- 
proach of Luc y & Sweenevl d!971l) indicates that there is a 66% 
probability that an eccentricity of or above the fitted value could 
have arisen by chance if the the underlying orbit is in fact circu- 
lar. As such, we impose a circular orbit, but we note that doing 
so has no signicant effect as the fitted eccentricity was so small. 

The median values and 1 cr uncertainties of the system pa- 
rameters derived from the MCMC model fit are presented in 
Table |4] The corresponding best-fitting transit light curves are 
shown in Figure [TJ and Figure [3] and the best-fitting RV curve is 
shown in Figure [2] 

Though the visual companion is 40 times fainter than WASP- 
3 1 and is resolved in the Euler and FTN images, we did correct 
the WASP photometry for the contamination prior to producing 
the MCMC solution presented. We checked the effect of the con- 
tamination by producing another MCMC solution using the non- 
corrected WASP photometry. The best-fitting parameter values 
were the same to within a tenth of an error bar. 

Without exquisite photometry, our implentation of MCMC 
tends to bias the impact parameter, and thus R„ and R v \, to 
higher values. This is because, with low-S/N photometry, the 
transit ingress and egress durations are uncertain, and symmet- 
ric uncertainties in those translate into asymmetric uncertain- 
ties in b and thus in R„. The effect on the stellar and planetary 
radii is larger for high-impact-parameter planets such as WASP- 
31b. Therefore we explored an MCMC with a main-sequence 
(MS) prior imposed (ICollier Cameron et al.l2 007). This employs 
a Bayesian penalty to ensure that, in accepted steps, the values of 
stellar radius are consistent with the values of stellar mass for a 
main-sequence star. The differences between the solutions with 
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Table 4. System parameters from the combined analysis 



Parameter (Unit) 


Value 


p CHI 

To (HJD) 

Tu (d) 

T n = Tta (d) 

a/R, 

b 

((") 


J.4UJ7U7U It U.UUUUUJ 

2455 192.6887 + 0.0003 
0.1103 ±0.0013 
0.0285 ± 0.0018 
8.00 ±0.19 
0.01615 ± 0.00027 
0.780 ±0.013 
84.41 ± 0.22 


K, (ms- 1 ) 

fl(AU) 

e 

TCORALIE (lUS 4 ) 

Atharps (ms _1 ) 


58.1 ±3.4 
0.04659 ± 0.00035 

(adopted) 
-124.923 ±0.035 
-5.268 ± 0.093 


M, (M ) 
R, (Re) 
P, (Po) 
logg, (cgs) 
T eS (K) 
[Fe/H] 


1.163 ±0.026 
1.252 ±0.033 
0.592 ± 0.042 
4.308 ± 0.020 

6302 ± 102 
-0.200 ± 0.090 


M v \ (M Jup ) 

^pl C^Jup) 

Ppl (pjup) 
ioggpi (cgs) 

(K) 


0.478 ± 0.029 
1.549 ±0.050 
0.129 ±0.014 
2.659 ± 0.036 
1575 ± 32 



and without MS priors are small and within errors, indicating 
that the transit light curves are of a quality such that the ingress 
and egress durations are measured sufficiently well. As such, we 
adopt the solution without the MS prior, which has slightly more 
conservative error bars. 
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Fig. 4. Colour-magnitude dia gram for WASP-3 1 a nd its compan- 
ion. Various isochrones from lMarigo et al. (2008) are given, with 
ages indicated in the figure. No de-reddening has been applied 
as A v = 0. 




5. System age 

Using 2MASS photometry we constructed a colour-magnitude 
diagram (Figure© for WASP-3 1 and its common proper motion, 
K-type companion. A distance modulus of 7.8 ±0.2 (360±30pc) 
is required to place the companion on the main-sequence, which 
puts WASP-3 1 between the zero-age main sequence and 1-Gyr 
age lines. Thus, if the two stars are of common origin, they 
are most probably around 1-Gyr old, with an approximate up- 
per age li mit of 4 Gy r. Combining the absolute V ma „ of an 
F7-8 star dGravl l2008h with the measured J - K and H - K 
colours of WASP- 3 1 from 2MASS and the intrinsic colours from 
iKoornneefi ([[983), we derive an interstellar extinction Ay = 
and a distance of 360 + 20 pc. This agreement with the dis- 
tance determined from the colour-magnitude diagram supports 
the inference that the companion is associated with WASP-3 1 
rather than being a mere line-of-sight neighbour. At a distance 
of 360 pc, the sky-projected separation of the two stars suggests 
that they would be separated by at least 12600 AU (0.2 light- 
year). 

Assuming aligned stellar-spin and planetary-orbit axes, the 
measured v sin / of WASP-3 1 and its derived stellar radius in- 
dicate a rotational period of P lot = 7.9 + 0.7 d. Comb ining this 
with the B — V colour of an an F8 star from iGravl (120081) . we used 
the relationship of lBarnesI d2007l) to estimate a gyrochronologi- 
cal age of 950 + 250 Myr. This is an upper limit as the star would 
be rotating faster, and so be younger, than suggested by spectro- 
scopic vsin/ if the stellar spin axis were inclined with r espect 
to the sky plane. We used the method of Max ted et alJ d201 ll) 



Fig. 5. Modified H-R diagram. The isochrones (Z = 0.012 w 
[Fe/H] = -0.20) for the ages 0.5, 1, 2, 3, 4, 5 and 10 Gyr and 
the evolutionar y mass tracks (Z _ .012 « [Fe/H] = -0.20; Y = 
0.30) are from lMarigo etalj (2008). To obtain the mass tracks, 
we performed a simple linear interpolation of their Z = 0.0008 
andZ = 0.017 tracks. 

to search for rotational modulation of the WASP light curves, 
as can be caused by the combination of magnetic activity and 
stellar rotation. No evidence of modulation was found. 

The lithium abundance (A T i = 2.75+0 1 0) found in WASP-3 1 
implies an age (ISestito & Ran dich 2005) between that of open 
clusters such as M34 (250 Myr; A Li = 2.92 + 0.13) and NGC 752 
(2 Gyr; Au = 2.65 + 0.13). However, lithium is a poor indicator 
of age for a star as hot as WASP-3 1 , and the measured abundance 
is consistent at the 1-cr level with that of the upper envelope of 
the 5-Gyr M67 (A Li = 2.55 + 0.18). 

We interpolated the stellar evolution tracks of Ma rigo et al.l 
(2008) using p, from the MCMC analysis and using T e g and 
[Fe/H] from the spectral analysis (Figure 0. This suggests an 
age of 4 + 1 Gyr and a mass of 1.0 ± 0.1 M Q for WASP-3 1 . 

Taking account of each age indicator, we suggest a likely age 
of 1^ 5 Gyr. 

6. Discussion 

With a mass of 0.48 M Jup and a radius of 1.55 fl Jup , WASP- 
31b has a density 13 per cent that of Jupiter and is ~0.3 Rj np 
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larger than predicted by standard models of irradiated gas giant s 
(iFortnev et al.ll2007h . Only WASP-17b (lAnderson et all 1201 Oh . 
which has a simil ar mass (0.49 Mj UD ), is known to have a lower 
density (0.06 pT„ r . lAnderson et al.ll201 lbl) . 

With an increasingly large sample of well-characterised 
planets, we can begin to make statistical infer ences as to 
the ph ysical reasons behind their diverse natures. Enoch et al. 
(1201 ll) showed the radii of 16 of the 18 known low-mass (0.1-0.6 
Mj llp ) planets strongly correlate with equil ibrium temperature 
and host-star metallicity. The calibration of lEnoch et alJ d201 ll) 
predicts a radius of 1.39 Rj np for WASP-31b. In a similar study, 
but using a different metallicity dependence and treating the 74 
known Jupiter-mass (0.2-2.5 Mj up ) planets, Anderson & Iro (in 
prep.) also found a strong correlation between planetary radius 
and equilibrium temperature and host-star metallicity. The cali- 
bration of Anderson & Iro (in prep.) predicts a radius of 1 .23 Rj up 
for WASP-31b. In each case, the predicted radius of WASP-31b 
is smaller than the measured radius (1.55 + 0.05 Rj U p). 

WASP-31 has a si milarly low metal licity to WASP-17 
([Fe/H] = -0.19 ± 0.09: lTriaude7ai1l2010l) . thus both WASP- 
31b an d WASP- 17b could reasonably be expect ed to have small 
cores dGuillot et al.ll2006t [Burrows et al.ll2007l) . However, this 
would only somewhat explain why the two planets are so large. 
Both planets are highly irradiated, with WASP- 17b being more 
irradiated than WASP-31b as, despite being in a slightly wider 
orbit (a = 0.052 AU), its host star is larger (R „ = 1.58 R Q ) 
and hotter (T eff = 6650 K; lAnderson et"aTll201 lbl) . This results 
in an equilibrium temperature for WASP- 17b hotter by 200 K 
than for WASP-31b and, from this, we could expect WASP- 
17b to be larger than WASP-31b. Both planets, though, are 
larger tha n predicted by stand ard models of irradiated giant plan- 
ets^. g. Fortn ev et al.ll2007l) . and by the empirical relations of 
lEnoch et alJ (1201 lb and Anderson & Iro (in prep.). Hence, it 
see ms likely that some addition al physics, such as Ohmic heat- 
ing (iBatvein & Stevensonll2010l) . is at play. 

The RV data place a stringent upper limit on WASP-3 lb's or- 
bital eccentricity (<? < 0.13; 3<x). It is therefore unlikely that tidal 
heati ng resulting from the cir cularisation of an eccentric orbit 
(e.g. iBodenheimer et ai]|2001l) was responsible for significantly 
inflating the planet. However, we could happen to be viewing 
the system soon after circularisation occurred and prior to the 
planet significantly contracting. This would have made finding 
the planet easier due to the greater transit depth. 

The metallicity of WASP-3 1 is at the lower end of what may 
be expected for a star of its age i n the Solar neighbou rhood (at a 
Galactocentric radius of 8.5 kpc; Magrini et al. 20091). 
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